Magnetic effects on nonlinear mechanical properties of a suspended carbon nanotube 
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We propose a microscopic model for a nanoelectromechanical system made by a radio-frequency 
driven suspended carbon nanotube (CNT) in the presence of an external magnetic field perpendic- 
ular to the current. As a main result, we show that, when the device is driven far from equilibrium, 
one can tune the CNT mechanical properties by varying the external magnetic field. Indeed, the 
magnetic field affects the CNT bending mode dynamics inducing an enhanced damping as well as 
a noise term due to the electronic phase fluctuations. The quality factor, as observed experimen- 
tally, exhibits a quadratic dependence on external magnetic field strength. Finally, CNT resonance 
frequencies as a function of gate voltage acquire, increasing the magnetic field strength, a peculiar 
dip-peak structure that should be experimentally observed. 



I. INTRODUCTION 

NanoElectroMechanical Systems (NEMS) made of 
suspended CNTs have received increasing attention 
recentl}M^. These devices are ideal for NEMS applica- 
tions due to the extreme mechanical properties of CNTs 
(low mass density and a high Young's modulus), resulting 
in a wide range of resonance frequencies for the funda- 
mental bending mode vibration (from MHz up to GHi^l 
range). In particular, CNT-based electromechanical de- 
vices working in the semiclassical regime (resonator fre- 
quencies in MHz range compared to an electronic hopping 
frequency from the leads of the order of tens of GHz) have 
attracted great interest due to t he extr emely large qual- 
ity factors {Q > 10^) attainabl c!^ ' ^ ' ^" '. Furthemore, at 
cryogenic temperatures, such devices behave as quantum 
dotsiSn^ with a strong interplay between single- electr on 
tunneling and bending mode mechanical motiorPE^HIll, 
This means that the electronic current is very sensitive 
to the CNT bending mode dynamics and can be used as 
a quantum (due to the intrinsic quantum nature of the 
charge carriers) measurement device. On the other hand, 
large quality factors allow to tune CNT mechanical prop- 
erties, e.g. the resonance frequency of the bending mode, 
by adjusting electronic parameters such as gate and bias 
voltages^'^^. As a main consequence, one expects that 
the application of a static magnetic field perpendicular 
to the CNT device modifies the above picture since the 
electronic current flow is affected by the Lorentz force. 
This, in turn, induces a change in the mechanical proper- 
ties of the CNT resonator that can be measured through 
the current itself. 

Motivated by recent experimentJ^Ell^ -^^e study a gen- 
eral model describing the mechanical properties of a 
radio-frequency driven suspended CNT-based NEMS in 
the presence of a transverse magnetic field. Indeed, we 
have recently shown that the effects of damping, spring 
stiffening and softening, and nonlinearity observ ed in 
similar devices in the absence of a magnetic 
can be successfully described in terms of a very sim- 



ple effective modeP^. In that model, t he CN T descrip- 
tion is reduced to a single electronic leveP^tSSl coupled to 
two metallic leads and interacting by means of a charge- 
displacement term with a single vibrational degree of 
freedom describing the bending mode. Due to the low 
frequency associated to the bending mod^^U^ the vibra- 
tional CNT dynamics is described by a classical Langevin 
equation^^ 22 28,^ ^j^jg simple approach, we were able 
to reproduce, in the absence of magnetic field, the main 
results reported in Ref.[2] and [3^ on a similar CNT- 
resonator device and even to predict features not yet 
observed. In particular, we have shown that when bias 
voltages are smaller than the broadening due to tunnel 
coupling, the resonance frequency shows a single dip as 
a function of gate voltage while, at bias voltages exceed- 
ing the broadening due to tunnel coupling, the resonance 
frequency shows a double dip structure. The successive 
experimental observation of the predicted effect h as 
confirmed the validity of the model adopted by n^^^ in 
describing these devices. 

The main result of the present paper is to prove that 
the external magnetic field modifies the mechanical prop- 
erties of the CNT-resonator. In particular, the magnetic 
field provides an additional damping mechanism for the 
resonator mechanical motion. Interestingly, a quadratic 
decrease of the quality factor Q as a function of the ex- 
ternal magnetic field strength, in quantitative agreement 
with the experiment performed in Ref.[T], emerges. 

When the device is driven far from the equilibrium 
we find that, increasing the magnetic field, the peculiar 
features (single or double dip) observed in the resonance 
frequency against gate voltage curves get distorted and 
acquire a dip-peak structure that could be experimentally 
observed. In particular, this effect should be more easily 
detectable at bias voltages that exceed the broadening 
due to tunnel coupling. 

We show that, at a fixed gate, bias voltage and temper- 
ature, if charge and current variations of the opposite sign 
occur, damping increases and Q is reduced. Vice-versa, 
charge and current variations of the same sign reduce 
damping with a consequent increase of quality factors. 
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FIG. 1. A carbon nanotube (CNT) subject to an external 
magnetic field H suspended between two normal metal leads 
biased by a voltage eVbias ■ 

As concerns the physical mechanism triggered by the 
magnetic field, we show that the application of a field 
perpendicular to the current flux modifies all the terms 
describing the CNT-resonator dynamics. Actually, the 
coupling with a transverse magnetic field introduces an 
electronic tunneling phase which depends on the me- 
chanical displacement of the CNT-resonator itself. This 
modifies the effective force acting on the resonator by 
a pure non-equilibrium correction term proportional to 
the magnetic field as well as to the electronic currenlPsl. 
Moreover, even at zero bias voltage, damping and dif- 
fusive terms are both modified by quantum electronic 
current-current fluctuations corrections whose strength 
is quadratic in magnetic field. Finally, we further show 
that, at zero bias, charge-displacement and magnetic field 
mediated electron-oscillator couplings cooperate behav- 
ing as a whole as a standard thermal bath at leads' tem- 
perature. 

The paper ends with a study of the device response 
when the CNT-resonator motion is actuated by an ex- 
ternal antenna at fixed frequency and amplitude. In 
this case, the current-gate voltage characteristic exhibits 
specific structures corresponding to the mechanical res- 
onances (antenna frequency equal to the bending mode 
frequency) . 

The paper is organized as follows: In Sec. II we dis- 
cuss the model. In Sec. Ill we construct, by means of the 
adiabatic approximation, the stochastic Langevin equa- 
tion for the dynamics of the oscillator including magnetic 
field and external antenna effects. In Sec. IV we present 
numerical results. 



II. MODEL 

We consider the system sketched in Fig. ([!]), which 
shows a single wall CNT suspended between two nor- 
mal metal leads. An external magnetic field H is applied 
perpendicular to the CNT. We also restrict the CNT me- 
chanical degrees of freedom to the low frequency bending 
mode and model it as a harmonic oscillator with fre- 
quency Wq. 

The electronic part of the CNT device is modeled as a 
single electronic level coup led to the leads through stan- 
dard tunneling termglSESI. The electronic Hamiltonian 



is 

(1) 

where CNT's electronic level has energy Ugate with cre- 
ation (annihilation) operators d^{d). The operators 
c|,^(cfc_a) create (annihilate) electrons with momentum 
k and energy ek,a — E^^a — /^a in the left {a = L) or 
right [a — R) free metallic leads. The chemical poten- 
tials in the leads /i^ and are assumed to be biased by 
an external voltage eVuas = P^l — Mfi- 

In the presence of a magnetic field, the phases of the 
tunneling ampl itude s with the leads depend on the CNT 
displacement a;,'22Ei] 

Vk"r{^) - Vk,Re'P\ (2) 

where p = deHL/2h gives the CNT momentum change 
induced by the Lorentz force when an electron tunnels 
from the CNT to a lead, and (5 ~ 1 is a numerical fac- 
tor determined by the spacial profile of the fundamental 
modd^. Above, e is the modulus of the electron charge, 
h is the Plank constant and L is the CNT length. Slow 
time varying tunneling amplitudes Vk{t) are also rele- 
vant in the c ase of adiabatic quantum pumping through 
quantum dot^^MS]^ Furthermore, when the external mag- 
netic field values are sufficiently small, such as the Zee- 
man splitting is negligible compared to broadening due 
to tunnel coupling, we can neglect the effect of the elec- 
tronic spin degrees of freedom (this issue will be consid- 
ered elsewhere^. 

For the sake of simplicity, we will suppose symmet- 
ric coupling Vk^L — Vk.R. E^nd a fiat density of states for 
the leads pk.a '-^ Pa, considered as thermostats at fi- 
nite temperature T, wi thin t he wide-band approximation 
iVk,a ^Va, a = L,Rf^. Definitely, the total tunnel- 
ing rate is hT = J2a=L,R^^a, with r„ = 2TTPa\Vk^a\''^ /h. 

The Hamiltonian of the mechanical degree of freedom 
is given by 

Hose = 7^ + \'m'Ujlx^i (3) 
Am I 

characterized by the frequency wq and the effective mass 
m {k — mu j^). The charge-displacement interaction is 
provided bji^sEsl 

Hint = A£n, (4) 

where A is the charge-displacement coupling strength and 
fi = d^d represents the electronic charge-density on the 
CN^. Definitely, the overall Hamiltonian is 

H = Hel + Hose + Hint- (5) 

For the experiment discussed in Ref.fl], one has a 
strong separation between vibrational (wq ~ 500 MHz~ 
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2/^eV) and electronic time scales (F ~ SOGhz) so that we 
can solve our model in the adiabatic limit, wq/F << 1. 
The experimental values of bias voltages and tempera- 
tures allow als o a sem i-classical treatment of the oscil- 
lator dynamic^^IIIZEIl, this paper, we will measure 
lengths in units of a;o = r, where r is a small frac- 
tion of CNT radius (r = 60pm) appropriate to resolve 
the CNT bending dynamics at relatively small temper- 
atures {T ~ 25mK). For the sake of simplicity, we 
will indicate dimensionless displacement variable with 
X. Energies are measured in units of hT — 200fieV, 
and times in units of Iq — I/loq. In terms of these 



units, the dimensionless spring constant is k/r 



1, 



since, following Ref.[T], the effective mass of the nan- 
otube is m = 1.3 x 10~'^^kg. Definitely, the adiabatic 
ratio is wq/F — 0.01, while the dimensionless temper- 
ature ksT = 0.01. Magnetic fields are measured in 
terms of the quantity B = jf- where our magnetic field 
unit is Hq = 2h/eLr ~ 16.6T, since the CNT length is 
L ~ 700nm. Throughout this paper, we keep fixed the 
dimensionless charge-displacement coupling strength to 
A = 0.1 (our force unit is hT/r), corresponding to an esti- 
mate Ep — A^/2fc ~ \^eV , implying a moderate coupling 
between the electronic and vibrational degrees of freedom 
(Ep/hujo = 0.5). Summarizing, the regime of the relevant 
parameters is hojo — Ep ~ fc^T << eV^/^{. < UT. 

In the next section, we show how adiabatic approxi- 
mation works on the coupled electron-oscillator problem 
in the presence of a transverse magnetic field. 



III. ADIABATIC APPROXIMATION 

As analyzed in the previous section, we work in the 
physical regime where the vibrational motion of the 
CNT-resonator is 'slow' with respect to all electronic en- 
ergy scales and can be considered "classical": ujq << F. 
This regime of the parameters leads to the adiabatic ap- 
proximation for the electronic problem. In contrast to 
previous works that treated the adiabatic approximation 
in the absence of a magnetic fielcP^^, we here investi- 
gate the effect of a transverse magnetic field on the elec- 
tronic problem described by the Hamiltonian Eq. ([5| . We 
remark that the adiabatic approximation has been used 
to describe larger system^ for the study of spectral and 
transport properties of organic semiconductor^2lMI]_ 



A. Adiabatic approximation for the electron 
problem in the presence of a magnetic field 

In this subsection, we show how the adiabatic approxi- 
mation on the electronic CNT level Green function works 
in the presence of a transverse magnetic field. 

Assuming a slow time dependence of electronic Green 
functions on the resonator displacement x, we are able to 
calculate truncated expressions for the CNT level Green 
functions which acquire a " slow ' time dependence and. 



at first order, a linear correction in the oscillator velocity. 
As a result of the adiabatic approximation, the truncated 
CNT level Green functions will depend on the instanta- 
neous value of the position and velocity of the resonator 
G''''''<'>(w,x,v). 

The adiabatic expansion of the Fourier transformed 
retarded CNT level Green function is 



G"'(w,a;,v) = G(o)(w,x) -I- G(i) (w, x, v), 
where the expression of G^g^(a;,x) is 



G[o,(w,a;) 



1 



and that of GJ^-^-) (w, a;, v) is 



(6) 



(7) 



G(i)(w,X, v) = lhUgate{x)G'^(o){uJ,x)- 



9G[o)(w,a;) 



dhu! 



(8) 



Above, Ugate{x) — Ugate + and the dot indicates the 
time derivative Ugate = A^ = Av. 

Using the adiabatic approximatioiP^ x{ti) — x{t2) — 
i(^o)(^i ~ ^2), we obtain for the lesser and greater com- 
ponents in Fourier space 

^feads (W, V) ~ ^f^ads,(0) ('^) + ^leads,(l) ^ v) (9) 

where the expression of '^f^^ds (o) ('^) 

Sws,(o)(^) = *[7irL/LM + nFfl/nH], (10) 

and that of ^ieads,{i)[<^,^) is 



■^leads,{l) 



w,v) = -leHi — 

y 0[eVMas\ 



(11) 

Above, we have defined H = 2ph/e. The adiabatic 
correction to the lesser component of the leads self- 
energy in Eq.(ll) is entirely due to the transverse mag- 



netic field and represent one of the main results of the 
present paper, in contrast to previous works where mag- 
netic field effects in t he adiabatic expansion were not 
considerecP^llIZlllIll. 

Definitely, for the CNT level occupation we get 

{n){x,v) ~ (n)(o)(a;) + (n)(i)(a;, v), (12) 
where at zero order in the adiabatic expansion we get 

dflLU 



ifL{io) + fR{Lo))C{uj,x), (13) 



with the spectral function C{uj^x) 
the CNT level given by 



2Q?G[o) (w, a;) of 



C{uj,x) 



{ri^~Ugate{xW + [flTY/A- 



(14) 
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The first order corrections in the adiabatic expansion are 
Unear in the oscillator velocity 



with 



(")(!) {X,v)= V[i?(i) (x) + i?(2) (x)] , 



R(2){x) 



r ./ 27r 
eH f dhuj 



g^^\u;)C{iu,x), 



(15) 

(16) 
(17) 



where we have defined the transmission function T{uj, x) 



T{u},x) 



and 



[{hoj-UgateixW + myAY 

d[fL{co) + fR{u:)] 



dfiijj 
d[eVbias] 



(18) 

(19) 
(20) 



Above, R{i){x) is the the adiabatic correction to the 
density related to the charge-displacement coupling de- 
scribed in the interaction Hamiltonian E q.([4l), alr eady de- 
scribed in many papers in the literaturePrnZm, R^^2){x) 
is the adiabatic correction exclusively due to magnetic- 
coupling effects modifying the electronic phase of elec- 
trons flowing from the leads to the CNT. 

Finally, in the hypothesis of symmetric coupling to the 
leads Tl = Yr, one can calculate the adiabatic expansion 
for the symmetrized current (/) = [{II) — (-^_r)]/2 

{i)M = l I ^|G'-(a.,a;)p(S«->(..,v)I]^-<(..,v) 



27r 

S^'>(w,v)I]«^<(a;,v)). 
Using Eqs.(|6]),([9]), we get 

(/)(a;,v)^(/)(o)(x)-h(/)(i)(x,v) 

where 

(/)(o)(.) = er/|^ 
with linear corrections in the oscillator velocity 

(/)(1)(X,V) =v[C/(i)(x) + C/(2)(x)], 

with 



(21) 
(22) 

{fd^)- fR{uj))C{u:,x), (23) 

(24) 



C/(i)(x) 



eA [ dfiLj 



g'^l\^)C{Lo,x)T{io,x), (25) 



Ui2){x) 
where we have defined 



(1),. _ djfU^) - fRjio)] ^ 

dflLO 

d[eVbtas\ 



(26) 

(27) 
(28) 



As already discussed referring to adiabatic corrections 
to the average charge density of the CNT level, C/(i)(a;) 
is the the adiabatic correction to the electronic current 
related to the charge-displacement coupling described in 
the interaction Hamiltonian Eq.Q, already described in 
Ref.[24] and [25]. U{2){x) is the adiabatic correction to 
the current exclusively due to magnetic-coupling effects. 
As we show below, Eq.([26| is one of the main result of 
the present paper, describing the increase of damping 
acting on the CNT resonator given by the application of 
a transverse magnetic field. 

In next subsection, we show that, even in the pres- 
ence of a transverse magnetic field, the dynamics of the 
CNT-resonator can be accurately described by a stochas- 
tic Langevin equation. 



B. Langevin equation for the oscillator 

In the absence of a magnetic field, the effect of the elec- 
tron bath and the charge-displacement coupling on the 
oscillator dynamics gives rise to a stochastic Langevin 
equation with a position dependent dissipation term and 
white noise forc^^. As in Ref. T5], even in the present 
case the equation for the oscillator dynamics can be writ- 
ten as follows 



mx + A{x)x = F(o)(a;) + D{x)$,{t) -|- A^xt cos{uJextt), 

(c(t)) = o, mm) = s{t-t'), (29) 

where ^{t) is a standard white noise term. We have in- 
cluded in our schematization the effect of an external 
antenna exciting the motion of the CNT, where A^xt, 
LOext represent the antenna amplitude and frequency, re- 
spectively. In this section, we describe how all the terms 
appearing in above equation modify in the presence of an 
external transverse magnetic field. 

The total force acting on the CNT-resonator is 



F 



-kx - X{n){x,Y) + H{i){x,Y). 



(30) 



The linear elastic force exerted on the oscillator is mod- 
ified by two relevant nonlinear correction terms: the 
former is proportional to the electronic charge density 
on the CNT level Eq.p^, while the latter to the elec- 



tronic current Eq. ( 23 1 . The first term, due to the charge- 



displacement interaction on CNT-resonator and propor- 
tional to A was already discussed in Refs.[2] and [15]. 
Far from equilibrium and in the presence of a mag- 
netic field, a magnetomotive coupling between the CNT- 
resonator displacement and the electronic flow through 
the device comes into play. Actually, the transverse mag- 
netic field introduces a phase in the electronic tunneling 
that is proportional to the displacement of the CNT res- 
onator as well as on the field strength. This originates 
a Lorentz-like additive correction, linear in the magnetic 
field strength and in the electronic current, to the average 
force acting on the resonator. 
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In the limit of the adiabatic approximation, the force 
Eq.( |30[ ) can be decomposed in different expansion terms. 
It expHcitly depends on the oscillator position x through 
Ugate [x] — Ugate + and velocity v. The force is 



F(x,v) =F(o)(2;) + F(i)(a;,v), 



(31) 



where 



F(o)(x) = ~kx - A(n)(o)(a;) + H{i)^o){x), (32) 

and 

F{i){x,^) = -A(n)(i)(a;,v) + jy(/)(i)(x, v) 

= -A{x)v. (33) 

The total damping term A{x) is given by three contribu- 
tions 



A(x)=A^{x)+A"{x)+A"'^{x), (34) 



where both 



A^(x) = Ai?(i)(x), 



(35) 



coming from the charge-displacement coupling, and 

A"{x) = -HU(^2){x), (36) 

due to magnetic field coupling, are positive definite. 
The function A^''^{x) is proportional to both charge- 
displacement coupling strength A and magnetic field H 



A"'\x) = XR^2){x)-HU^,){x), 



(37) 



and is not positive definite. Remarkably, we have verified 
that the whole sum appearing in Eq. ( 34 ) is positive def- 



inite in all parameters regime of our model. This shows 
that, using a spinless fermionic model in the presence of 
normal (not ferromagnetic) electronic leads, the CNT- 
resonator experiences no negative damping regions. This 
is in contrast to results of Ref. ^3.1j, where the authors 
use a normal and a ferromagnetic lead and observe neg- 
ative damping and consequent nano-electromechanical 
self-excitations of the CNT-resonator system. 

A fluctuating term has to be be included to take cor- 
rectly into account the effect of the bath degrees of free- 
dom. When a magnetic field is present, the force-force 
fluctuations are given by three contributions (see ap- 
pendix A) 



{SF{t)SF{t')) = \^{5h{t)5n{t')) + 
- H\[{Sh{t)Si{t')) 
+ H^{5i{t)5i(t')), 



{6i{t)5n{t'))] 



(38) 



where we get a mixed current-density fluctuation con- 
tribution [{5n{t)5i{t')) + {5i{t)6n{t'))], and a current- 
current fluctuation contribution {SIa{t)6I{t')) to the 



In the adiabatic limit, exploiting the effect of the 'fast' 
electronic environment on the oscillator motion, one de- 
rives 



{SF{t)SF{t')) = D{x)6{t-t'), 
where in the presence of a magnetic fleld we have 
D{x) = D^{x) + D"{x) + D^'^(x), 

with 

D\x) = X'h I ^G<)(c.,x)G>)(c.,x) = 



(39) 



(40) 



(41) 



and 



D"\x) = fx 



x(I]f,f(w)s[„<(^) 



eXH 



did 



^\Gl,^iu;,x)\^Ciu;,x)x 



27r [{hLd-UgateixW + myA]^ 

x{[hTLrfL{u;){l - /i(w)) - [nrfl]Vfl(w)(l - fRiuj))} 

(42) 

where C{uj,x) — — 2SyGL> (w, a;) is the electronic spectral 



function of the electronic level defined in Eq.(14|. The 



noise strength contribution coming from current-current 
fluctuations is 



{l~T{uj,x)) + [fL{^){l- fL{^)) 



+ 



(43) 



where r(a;, x) 



C{uj,x). In the absence of electron 



bias voltage, one has D{x) = 2kBTA{x), that is the 
fluctuation-dissipation condition is verifled for each flxed 
position X. Moreover, it is possible to show that in the 
chosen units, the dimensionless damping A{x) (Eq.(34)) 
and diffusive term D{x) (Eq.(40)) result proportional to 
the adiabatic ratio ujq/T. 

It is important to point out that, when there is no 
intrinsic charge-displacement coupling (A = 0), in the 
absence of the antenna effects and at zero bias {Vbias = 
0), the oscillator is still governed by a Langevin equation 

mx -I- A"{x)x ^kx + ^D"{x)^{t), (44) 

with a harmonic force F(^q^(x) — —kx, an intrinsic 
positive-definite dissipative term A^{x), and a diffusive 
term (x) proportional to the thermal current-current 
noise. Looking at Eqs.(36) and (43), one can clearly see 
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FIG. 2. (Color online) Panel(a): minimum of the effective 



potential (coming from the force Eq.(32l) affecting the CNT- 
resonator as a function of the magnetic field at low bias 
eVbias = O.lfir {Vbias ~ in our units). Solid (black) 

line indicates Ugate = 0, dashed (red) line Ugate = 0.45 
{Vbias ~ 90/1 V in our units), dotted (green) line Ugate = 1.0 
{Vbias = 200^1^ in our units). Panel(b): same as above at 
large bias eVbias = l.SftF {Vuas = SOO/iV in our units). Solid 
(black) line indicates Ugate = 0, dashed (red) line Ugate = 0.75 
(Vbias = 150^F in our units), dotted (green) line Ugate = 1.5 
{Vbias = 300^1 in our units). 



that a natural quadratic dependence of damping and dif- 
fusive strength on the magnetic field emerges. This can 
be explained observing that, even at zero bias voltage, 
the electronic tunnehng events, whose phase is depen- 
dent linearly on the CNT displacements as well as on 
the magnetic field strength, perturb the CNT mechani- 
cal motion with a force with zero average (due to (/) = 0, 
H can be also different from zero) and square mean pro- 
portional to the magnetic field square. Definitely, even 
in the absence of external bias voltage Vbias, the mag- 
netic field applied perpendicular to the CNT couples to 
the bending mode dynamics behaving as a surrounding 
thermal bath at leads temperature ksT . 

We end this section with a systematic study of the 
spatial dependence of the total force, the damping (see 
Fig. ([3])) and diffusive terms (see Fig. ([4])) as a function of 
the bias voltage as well as on the magnetic field. 

As concerns the total force acting on the CNT res- 
onator, we point out that, for the magnetic field strengths 
investigated in this paper, the effective potential pre- 
serves its parabolic shape with a displaced minimum and 
renormalized curvature. For instance, when a left-to- 
right current fiows through the device (see the sketch 
in Fig.([T])) in the presence of a positive magnetic field 
(outgoing from the sketch reported in Fig.([T])), the CNT- 



FIG. 3. (Color online) Spatial dependence of the dimension- 
less damping coefficient A{x) at low bias (Panels(a-b-c)) and 
at large bias voltage applied (Panels (d-e-f)). See main text 
for discussion. 



resonator effective potential minimum is displaced to- 
wards positive displacements x with respect to the min- 
imum set by the charge-displacement interaction (see 
Panel (a-b) of Fig.([2])). In Fig.([2]), one can observe that, 
the minimum of the effective potential acting on the res- 
onator depends linearly on the magnetic field strength. 
This comes from the linear dependence on the mag- 
netic field of the Lorentz-like correction term to the force 
Eq.(|32|. In particular, as shown in Panel(a) in Fig.(|2|, 
in the low bias regime for the device (small compared 
to the broadening of the CNT level), the larger is the 
gate voltage, the smaller is the displacement of the po- 
tential minimum as a function of the external magnetic 
field with respect to the shift produced by the charge- 
displacement interaction on the CNT (whose position is 
indicated by a (black) square for Ugate = 0, a (red) cir- 
cle for Ugate = 0.45 {Vbias = 90/^^ in our units), and 
a (green) triangle Ugate = 1.0 {Vbias = 200/xl/ in our 
units)). This can be explained observing that in the low 
conducting regime of the device the resonator is less effec- 
tively coupled with the electronic subsystem. In the large 
bias regime (Panel(b) of Fig. ([2])), a smaller magnetic field 
is sufficient to displace the potential minimum of the 
same quantity produced by the sole charge-displacement 
interaction on the CNT. Again, the larger is the gate 
voltage, the smaller is the displacement of the potential 
minimum as a function of the external magnetic field with 
respect to the shift produced by the charge-displacement 
interaction on the CNT. 

The renormalization of the effective potential curva- 
ture, that is of the resonance frequency of the resonator, 
will be discussed in subsection B of next section. 
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In this section, we limit ourself to discuss the damp- 
ing term A(x), since for the diffusive term D(x), un- 
less explicitly stated, a similar analysis can be done. As 
shown above (see Eq.(34)), we can distinguish between 
three contributions to the friction affected by the oscil- 
lator: a pure charge-displacement contribution A'^{x), 
depicted in Fig.(|3]) with a solid (black) line, already 
discussed in Refs!^2j and [55]; a damping contribution 
due to current-current fluctuations A^{x), depicted in 
Fig. ([3]) with a dotted (blue) line; a mixed damping term 
due to current-density fluctuations (not positive defi- 
nite), indicated by A-'^^^x) and depicted in Fig.([3j) with 
a short-dashed (red) line. The total damping A{x) is re- 
ported with a dashed (pink) line. As one can observe 
in Panel(a), at low bias voltage, when the external mag- 
netic field strength is smaller than charge-displacement 
coupling A, the damping contributions coming from the 
current-current A^ {x) and current-density A^^ {x) fluc- 
tuations are negligible with respect to that generated 
by the pure charge-displacement contribution A^{x). In 
Panel (a), A'^{x) and A^ (x) have a single peak struc- 
ture centered at a; — 2Ugate/^ — 0, while A^^ {x) is an 
odd symmetric function with respect to this point. We 
point out that these peculiar structures emerge only at 
larger values of magnetic field (Panels(b-c) of Fig.(|3|). 
The total damping affecting the resonator is peaked at 
configurations where large density variations take place 
\x - 2Ugate/M < TtA:/\. Indeed, the density of the CNT 
level goes from a region x — 2Ugate/^ < —TiT/\ corre- 
sponding to almost completely filled states ((n) ~ 1) to 
a region x — 2Ugate/^ > hT/X corresponding to com- 
pletely empty states {{h) ~ 0). Definitely, the CNT 
level experiences an unit c harge variation across the 
\x - 2Ugate/M < 7ir/A regioij™. 

At large bias voltages applied (Panel (b)), A^{x) has 
two peaks centered at a; — 2Ugate/X ~ eVbias/2A and 
X — 2Ugate/^ — — eVf,ias/2A, respectively. A^ {x) shows 
the same behavior, while A^^ {x) is an odd symmet- 
ric function with respect to these two points. As in 
Panel (a), A^{x) and A^^ [x) are negligible with re- 
spect to A^{x). The total damping affecting the res- 
onator is peaked at configurations where the CNT level 
experiences a half- unit charge variat ion ac ross the \x — 
2Ugatel\ ± eVuasl2\\ < fiT / X region^™. 

When the external magnetic field is turned on, an en- 
hanced damping as well as as noise strength emerges with 
a quadratic dependence on the magnetic field intensity 
Eqs.([34])-(|40|. In Panels(b-e) of Fig.([3|, one can observe 
that, as the dimensionless ratio B/A is equal to one, the 
total damping affecting the resonator is only slightly per- 
turbed by the application of the magnetic field. At low 
bias, A{x) preserves its single peak structure with an 
enhanced strength (dashed (pink) curve in Panel (b) of 
Fig.([3|). At large bias, the strength of the two peaks 
becomes asymmetric, with an enhanced damping of the 
peak sX X — 2Ugatel^ — ~eVhiasl2. This effect can be 
explained as follows: when a magnetic field is applied to 
the device, the resonator starts to feel even the varia- 
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FIG. 4. (Color online) Spatial dependence of the dimension- 
less diffusive coefficient D{x) at large bias voltage applied 
(Panefs(a-b-c)). See main text for discussion. 



tions of the electronic current flowing through the CNT 
as a function of the gate voltage (see Eq.(32)). These 



current variations are positive for x — 2UgateTX < and 
negative otherwise. At x ~ 2Ugate/X ~ ~eVbias/2, large 
negative variations of the electronic density and positive 
variations of the electronic current cooperate giving an 
enhanced damping. 

We intend now to study the regime realized when the 
external magnetic fleld strength is larger than charge- 
displacement coupling strength A. In this case, the 
contribution to the damping coming from the current- 
current fluctuations A^ (x) are dominant with respect to 
those corresponding to charge-displacement A^{x) and 
current-density A^^ (x) fluctuations. In the low bias 
regime, the total damping term preserves its single peak 
structure which, due to the intrinsic asymmetry of the 



current-density term A 



\H( 



is slightly distorted. For 



the same reason, in the large bias regime, the double dip 
structure of the total damping term is preserved with 
an enhanced asymmetry. In the large magnetic fleld 
regime, it is important to point out the particular spatial 
dependence of the noise strength D{x) (Panels (c-f) of 
Fig.Q). Here, the noise contribution due to the current- 
current fluctuations {D^ (x)) emerges with the charac- 
teristic double peak structure even at low bias regime 
(dotted (blue) curve in Panel (c) of Fig. Q ). Comparing 
the dashed (pink) curves in Panel (c) of Figs.(|3])-(|4|, one 
can observe that, even at low bias voltage, the applica- 
tion of a large magnetic field drive the CNT-resonator far 
out of equilibrium, breaking the validity of the Einstein 
relation D{x) — 2kBTef fA{x) with an effective tempera- 
ture. Far from equilibrium, this relation is strictly valid 
only at very low bias voltage o^^'^^ ' . 
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In the next section, we study numerical results of 
our model concerning mechanical properties of CNT- 
resonator (resonance frequency and quality factor) as well 
as the electronic observables inherent to the transport 
problem (I-V characteristic). 



IV. MECHANICAL AND ELECTRONIC 
CHARACTERISTICS OF THE DEVICE 

Given the assumption about the separation between 
the slow vibrational and fast electronic (tunneling) 
timescales, the problem of evaluating a generic observ- 
able (electronic or not) of the system reduces to the eval- 
uation of that quantity for a fixed position x and velocity 
v of the oscillator, with the consequent averaging over the 
stationary probability distribution P(a;, v). From the so- 
lution of the Langevin equation ( 29 1 , one can determine 



the distribution P{x, v) which allows to calculate all the 
electronic observables O: 



(O) 



dxdvP(x, v)0{x, v). 



(45) 



We analyze in the next section the effects of the magnetic 
field on the mechanical as well as electronic properties of 
the device. 
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FIG. 5. (Color online) Device quality factor as a function the 
gate voltage Ugate for different magnetic field strengths at 
low (Panel (a)) and large (Panel (b)) bias voltage. Panels(c- 
d) Same as above for the average total damping {A{x)) of the 
system. Panel(a-c): solid thick line B = 0.0, solid normal- 
thickness line B — 1.5, and solid thin line B = 3.0. Panel(b- 
d): solid thick line B = 0.0, solid normal-thickness line B = 
0.2, and solid thin line B = 0.4. 



A. Device Quality factors 

One of the main findings of Ref.[T] is the observation 
of a quadratic dependence of the device quality factor Q 
on external magnetic field strength. Within our model, 
as also stressed in the previous sections, such a quadratic 
dependence on B emerges naturally. In order to include 
back-actions effects of the out of equilibrium electronic 
bath on the resonator, we have calculated the average 
device quality factor as 

Q= / dx^P{x), (46) 

where A{x) is the total damping at a particular resonator 
displacement x and P{x) is the reduced displacement dis- 
tribution probability of the CNT-resonator. We have ver- 
ified that this particular way of extracting quality fac- 
tors is completely equivalent to measure the width at 
half-high in the current-frequency curves obtained in the 
linear response to an external antenna exciting the nan- 
otube motiorfi^. 

Motivated by the experiment performed in Ref . [1] and 
by recent experimental study on a similar CNT devictPl, 
we here performed a systematic study of the quality cal- 
culated from our model as a function of the bias, gate 
voltage as well as on the magnetic field. In Panels (a- 
b) of Fig. ([5]), we investigate the device quality factor Q 
as a function of gate voltage in the low and large bias 



voltage regime, respectively. In the absence of a trans- 
verse magnetic field, we reproduce the qualitative behav- 
ior obtained in the experiment of Ref.fl^. When bias 
voltages are smaller than the broadening due to tunnel 
coupling, the quality factor shows a single dip feature 
(solid (black) thick line in Panel (a) of Fig.([5])). At bias 
voltages that exceed (or are equal to) the broadening 
due to tunnel coupling, the quality factor shows a dou- 
ble dip structure (solid (black) thick line in Panel (b) of 
Fig. ([5])). This behavior, already addressed in Refs. [T4] 
and jl5 , can be easily explained looking at the average 
charge and dissipation of the CNT-resonator. As also 
discussed referring to total damping affecting the CNT- 
resonator in the previous section, at low bias voltage and 
in the absence of magnetic field, the total average damp- 
ing affecting the resonator is peaked at electronic config- 
urations where the CNT level experiences an unit charge 
variation across the region where the small conduction 
window is placed |J7gQte| < (solid (black) thick line in 
Panel (c) of Fig. ([s])). At large bias voltages, the conduc- 
tion window, whose extension is proportional to eVbias, 
becomes larger than the broadening of the CNT level, so 
that the total average damping affecting the resonator 
is peaked at electronic configurations where the CNT 
level experiences half-unit charge variations, that is at 

\Ugate - eVMas/2\ < HT and \Ugate + eVbiasf^l < hT. 

When the transverse magnetic field is turned on, the 
above scenario modifies as follows. At low bias volt- 
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ages, the total damping affected by the CNT-resonator 
increases quadratically with the field at every point in the 
configuration space of the oscillator. Moreover, the CNT- 
resonator distribution probabilities P{x) depend slightly 
on the magnetic field as well as on the gate voltages and 
are actually centered at configurations close to the har- 
monic potential minimum a; ~ in the absence of charge- 
displacement interaction A. The overall result is an en- 
hanced average total damping as one increases the mag- 
netic field (solid normal-thickness (J5 = 1.5) and thin 
{B = 3.0) (black) lines in Panel (c) of Fig.([5])) and a 
corresponding decrease of the quality factor in all the 
gate voltage range investigated (solid normal-thickness 
{B = 1.5) and thin {B = 3.0) (black) fines in Panel 
(a) of Fig. (5])). At large bias voltage, the P{x) still 
depends only slightly on the magnetic field but is very 
spread on the configuration space. Therefore, the aver- 
age in Eq.(46) reproduces the spatial dependence struc- 
ture of the total damping coefficient reciprocal 1/A{x). 
The double peak structure of the average total damping 
term (solid thick (black) fine in Panel (d) of Fig. ([5])) is 
canceled by the magnetic field, giving a single peak at 
Ugate = eVbias/'i whcrc a cooperation between negative 
charge-density and positive current variations take place 
(solid normal-thickness {B — 0.2) and thin {B — 0.4) 
(black) lines in Panel (d) of Fig.([5|). As a consequence, 
the quality factor loses its double dip structure getting a 
single dip at Ugate — —eVhiasj'^ (sofid normal-thickness 
(B = 0.2) and thin (B = 0.4) (black) lines in Panel (b) 
of Fig.(l5l) 



We intend now to study the device quality factors as 
a function of the transverse magnetic field B comparing 
different conducting states of the device. In Panel (a) of 
Fig. ([6]), one can observe calculated device quafity factors 
as a function of the magnetic field at a low conducting 
state of the device {Ugate = 0.45). Different curves, from 
thicker to thinner, refer to increasing bias voltages ap- 
plied eVbias = 0.1 — 0.75 — 1.5. At every fixed bias volt- 
age, a clear quadratic dependence of the total average 
damping on the magnetic field strength is observed (not 
shown in Fig.(ro|), with a Lorentzian shape of the quality 
factor curves (see Panel (a) of Fig.Q). It is important 
to point out that the range of magnetic field strengths 
experimentally investigated in Ref . llj , = — 3T, cor- 
responds to small magnetic fields in our units (wc recall 
that Hq = 16.6T). Remarkably, at low bias and small 
magnetic fields, a quadratic decrease of the Q against 
magnetic field is observed (see solid thick (black) line in 
Panel (a) of Fig.(|6|). In Fig. ([7]), we show the quantitative 
agreement between experimental and calculated quality 
factors against magnetic field when the device is in a low 
conducting state, with eVbias — 0.1 and Ugate = —0.45. 
The slight increase of the quality factor Q as a function 
of the field for small magnetic fields, is due to asymmetry 
introduced by the gate voltage Ugate = 0.45 applied to 
the device (see also Panel (a) of Fig. ([6])). For gate voltage 
equal to zero, that in the high conducting state of the de- 
vice, the calculated Q against B curve is a parabola with 
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FIG. 6. (Color online) Device quality factor as a func- 
tion tlie magnetic field strength B for different bias volt- 
ages (Solid thick line eVuas = 0.1, solid normal-thickness 
line eVbias = 0.75, and solid thin line eVbias ~ 1-5) at low 
(Panel (a)) and high (Panel (b)) conducting states. Panels(c- 
d) Device quality factor as a function the bias voltage eVbias 
for different magnetic field strengths at low (Panel (c)) and 
high (Panel (d)) conducting states. Panel(c): solid thick line 
B — 0.0, solid normal-thickness line B = 0.05, and solid thin 
line B = 0.1. Panel(d): solid thick line B = 0.0, solid normal- 
thickness line B — 0.25, and solid thin line B — 0.5. 
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FIG. 7. (Color online) Device quality factor as a function 
the magnetic field strength. Squares represent experimental 
values obtained in Ref. 1 at T = 25mA', Vbias = O.SmV 
and distance from the current peak Ugate = —90ml/. Solid 
(red) line is calculated Q at ksT = 0.01, eVbias = 0.1 and 
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0.45. 
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a maximum at zero magnetic field applied. 

Coming back to Panel (a) of Fig. ^ , one can observe an 
interesting increase of the quality factor peak as a func- 
tion of the bias voltage. In particular, for eVuas = 1.5 
(thinner line in Panel (a) of Fig. ([6])) a quality factor peak 
at i? ~ Ugate/"^ — 0.225 occurs. This can be directly re- 
lated to the average total damping dip, not shown in 
Fig.Q. This effect can be explained noting that, when 
the bias voltage applied to the electronic device is in- 
creased, a transition from a single peak to a double peak 
structure in the spatial dependence total damping af- 
fected by the CNT-resonator con be observed (compare 
Panel(a) and (d) of Fig.(|3])), while at the same time, the 
displacement distribution probabilities P{x) spread on 
the configuration space remaining centered at configura- 
tions close to the harmonic potential minimum a; ~ 
characteristic of the low bias regime. The overall result 
is a reduction of the average total damping affecting the 
CNT-resonator whose minimum is translated by a quan- 
tity proportional to the gate voltage applied to the de- 
vice. This argument becomes even more clear when no 
gate voltage is applied to the device which is therefore 
placed in a high conducting state. In this case a perfect 
symmetry of Q-factor curves with respect to zero mag- 
netic field is obtained (see Panel (b) of Fig.Q). 

We end this section with a study of the device qual- 
ity factors as function of the bias voltages and magnetic 
fields comparing low and high conducting states of the de- 
vice. In Panel (c) of Fig.([6|, one can observe calculated 
device quality factors as a function of the bias voltages at 
a low conducting state of the device {Ugate = 0.45). Dif- 
ferent curves, from thicker to thinner, refer to increasing 
magnetic field applied to the device B = 0.0 — 0.05 — 0.1. 
At zero magnetic field, a clear double dip feature in the 
quality factor Q, as experimentally observed in Ref. |14j . 
is visible. This can be explained looking at the aver- 
age total damping and in terms of the average charge 
present on the CNT level. The total average damp- 
ing, in the absence of magnetic field, has two peaks at 

eVbias = -^Ugate = "0.9 and at eVbias = "^Ugate = 0.9. 

Indeed, as also discussed previously, the total average 
damping is peaked at electronic configurations where the 
CNT level experiences half-unit charge variations, that 
is at \Ugate- e.Vuas/'^\ < ?ir and lUgate + eVbiasf^l < ?ir. 
Therefore, when the edges of the conduction window 
(whose width is proportional to eVbias) meet the CNT 
level energy (given by Ugate), a maximum total average 
damping (minimum quality factor) is observed. As also 
discussed in reference of Fig. ([5]), the double peak struc- 
ture of the average total damping term is canceled by the 
magnetic field, giving a single peak at eVbias = —^Ugate 
where a cooperation between negative charge-density and 
positive current variations take place. As a consequence, 
the quality factor loses its double dip structure getting 
a single dip at eVbias = '^'2Ugate (solid normal-thickness 
{B = 0.05) and thin {B = 0.1) (black) lines in Panel (c) 
of Fig. ([6])). 

In Panel (d) of Fig. ([6]), we show calculated device qual- 




FIG. 8. (Color online) Panel (a): Resonator frequency against 
eflective gate voltage calculated as minimum of the effective 
potential in the static approximation at small bias eVbias = 
0.1 for different magnetic field values: Solid thick line B = 
0.0, solid normal-thickness line B — 1.5, solid thin line B = 
3.0. Panel (c): same as Panel (a) at large bias eVbias = 1.5?ir 
for different magnetic field values: From thicker to thinner 
line B = 0.0 - 0.1 - 0.2 - 0.3 - 0.4. Panel (b-d): Resonator 
frequency calculated using an external antenna (with A^^t = 
10"'') at mechanical resonance against effective gate voltage 
for same parameters of Panel (a-c), respectively. Dashed (red 
online) and solid (blue online) portions of each curve indicate 
resonance frequency values with positive and negative current 
change A/, respectively. In Panel (d), only curves referring 
to magnetic field strengths B = 0.0 — 0.2 — 0.4 are reported. 



ity factors as a function of the bias voltages at a high 
conducting state of the device {Ugate = 0.0). Different 
curves, from thicker to thinner, refer to increasing mag- 
netic field applied to the device B = 0.0 — 0.25 — 0.5. 
As above, at zero magnetic field, a single dip feature 
in the quality factor Q, as experimentally observed in 
Ref. |14j ■ is visible. This behavior can be discussed with 
the same argument given for discussing the Panels(a-b) 
of Fig.(|6|, where a reduction of the total average damp- 
ing as a function of the bias voltage applied to the device 
was observed. Again, a decrease of the quality factor in 
all the gate voltage range investigated as a function of 
the magnetic field is observed (see Panel (d) of Fig.Q). 



B. Resonance frequency renormalization and 
current-voltage curves 

In this section, we address the magnetic field effects on 
the renormalization of the CNT-resonator resonance fre- 
quencies and its back-action effects on the current volt- 
ages curves of the device. In order to study the CNT 
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resonance frequency renormalization as function of the 
gate voltage, we have compared results coming from two 
ways of evaluation of the resonance frequencies. In the 
first method, referred to as Static, we evaluate the po- 
sition of the minima of the static potential arising from 
the generalized force acting on the resonator (Eq.([32|). 
In the second method, referred to as Dynamic + antenna, 
we have analyzed, at every fixed value of the gate volt- 
age, all the traces of electronic current as a function of 
the antenna frequency reporting with a red (blue) dot 
the resonance frequency values with positive (negative) 
current change AI = I — Iq with respect to background 
value Iq obtained in the absence of the antenna. 

In Fig.Q, we report the resonance frequencies of the 
CNT-resonator as a function of the gate voltage compar- 
ing the two methods outlined above. We address the low 
bias regime in Panels (a-b), while the large bias regime is 
investigated in Panels(c-d). In Panel (a) of Fig.(|8]), dif- 
ferent curves, from thicker to thinner, refer to increasing 
magnetic field applied to the device B = 0.0 — 1.5 — 3.0. 
The same description was done in Panel (c), where differ- 
ent curves refer to increasing magnetic field in the range 
B = 0.0- 0.2 - 0.4. The thicker (black) lines in Pan- 
els (a) and (c) , corresponding to the absence of magnetic 
field, reproduce qualitatively all results experimentally 
observed in Ref . [T3] : when bias voltages are smaller than 
the broadening due to tunnel coupling (Panel (a)), the 
resonance frequency shows a single dip as a function of 
gate voltage. At bias voltages that exceed the broad- 
ening due to tunnel coupling (Panel (c)), the resonance 
frequency shows a double dip structure. Actually, in this 
regime, the onset of a double dip structure was already 
predicted by us in Ref. (15]. It is important to point out 
that the resonance frequency renormalization curves ob- 
tained in the presence of the external antenna (Panel 
(b-d) of Fig.Q) have the same qualitative behavior (as 
a function of the gate) of those obtained in the static 
approach. In the presence of an external antenna with 
a finite amplitude, renormalization effects in the reso- 
nance frequencies are less pronounced due to nonlinear 
softenin^2li5]_ 

As already analyzed in Refs.[5] and (TS], when the de- 
vice is in a low current-carrying state, a peak in the 
current-frequency curve signals the mechanical resonance 
(whose position is indicated by thin (red) lines in Panels 
(b-d) of Fig.([8])), while in a high current-carrying state, 
a dip in the current-frequency curves is observed (whose 
position is indicated by thick (blue) lines in Panels (b- 
d) of Fig.Q). In the presence of a transverse magnetic 
field, the different character of low and high conducting 
states, signaled by a peak or a dip in current-frequency 
curves is preserved (curves (2) and (3) in Panels (b-d) of 
Fig.®). 

The peculiar features of CNT-resonator frequency 
renormalization as a function of the gate can be ex- 
plained with the same argument used to describe the 
quality factors behavior in the previous section. Indeed, 
in the absence of magnetic field, the resonator frequency 




FIG. 9. (Color online) Panel (a) Average electronic current 
flowing through the CNT level at low bias {eVbias = 0.1) as 
function of the gate voltage for different values of the mag- 
netic field and in the presence of a external antenna applied 
to the device at fixed frequency LJext = 0.9975 and amplitude: 



solid (black) line 



10"^, (dashed (red) line B = 1.0 



Aext = 10"", dotted (green) line B = 1.5 Ae^t = 10"^ '^, 
dashed-dotted (blue) Une B = 1.0 A^^t = lO"^'^). Panel (b) 
Average electronic density on the CNT level for the same pa- 
rameter values as in Panel(a). See the main text for detailed 
discussion. 



renormalization is maximum at electronic configurations 
where the the CNT level experiences the largest charge- 
density variations against the gate voltage. 



k + X 



d{f 



dU, 



gate 



(47) 



Actually, at low bias voltages, an unit charge density 
variation across the region where the small conduction 
window is placed iC/natel < KT (solid (black) thick line 
in Panel (a) of Fig. ([8])) occurs. At large bias voltages, 
the CNT frequency renormalization is larger at electronic 
configurations where the CNT level experiences a half- 
unit charge variation, that is at \Ugate — s.Vuas/'A < ^-^ 
and at \Ugate + eVbias /2| < fiT. 

When the transverse magnetic field is turned on, the 
above scenario modifies as follows. The resonator fre- 
quency renormalization is larger at electronic configura- 
tions where the CNT level experiences the largest charge- 
density and current variations against the gate voltage. 
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gate 
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gate 



(48) 



At low bias voltage, the single dip feature in the 
CNT-resonator resonance frequency gets distorted (solid 
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normal-thickness line in Panel(a) of Fig.Q) and ac- 
quires, in the limit of large magnetic field (solid thin 
line in Panel(a) of Fig.Q), a dip-peak structure that 
could be experimentally observed. Actually, the peak 
observed at Ugate — 0.3 corresponds to an hardening of 
the CNT-resonator resonance frequency. This effect can 
be explained as follows: when a magnetic field is applied 
to the device, the resonator starts to feel even the varia- 
tions of the electronic current flowing through the CNT 
as a function of the gate voltage (see Eq.(48)). These 
current variations are positive for Ugate < and neg- 
ative otherwise. At Ugate — 0.3, the positive (due to 
the positive sign of the magnetic field) variations of the 
electronic current overcome the negative variation of the 
electronic density giving a hardening in the CNT reso- 
nance frequency. At Ugate — —0.25, one has negative 
variations of both density and current, obtaining a more 
pronounced softening in the resonance frequency. The ef- 
fect outlined above is more pronounced in the large bias 
regime (see Panel (c) of Fig.(|8|). Here, the magnetic field 



gives an enhanced softening dip at Ugate — —eVbias/'^ 
and an hardening peak at Ugate — s^UasI^^ where pos- 
itive variations of the electronic current cooperate with 
negative variation of the electronic density. In both low 
and high bias regime, the hardening effect outlined above 
could be experimentally observed. 

The peculiar renormalization frequency effects dis- 
cussed above have a nontrivial back-action effects on the 
electronic density and current-gate voltage characteristic 
of the device. In Fig.Q, we study the electronic CNT 
level density and current as a function of the gate volt- 
age in the presence of an external antenna at fixed am- 
plitude and frequency ixiext = 0.9975 (corresponding to 
the horizontal line in Panel (b) of Fig.(|8])) in the low 
bias regime of the devic^^. When the external antenna 
frequency becomes equal to the proper frequency of the 
resonator, we observe a dip structure in both density and 
current at a gate voltage corresponding to high conduct- 
ing states of the device (solid (black) line in Panels (a-b) 
of Fig.Q). This feature, that could be experimentally 
observed, is considered as a "dip" with respect to cor- 
responding curves in the absence of antenna or with an 
antenna frequency far from the range of the proper fre- 
quencies of the CNT-resonator (not shown in Fig.(|9|). 

When a transverse magnetic field is applied to the 
device, the CNT frequency renormalization profile as a 
function of the gate voltage changes (see Fig.Q). There- 
fore, the mechanical resonance condition between the ex- 
ternal antenna frequency and the proper frequency of the 
resonator occurs at different electronic gate voltages. For 
sufficiently large magnetic fields, the resonance can occur 
in correspondence of a low conducting state of the de- 
vice. As one can observe in Panels (a-b) of Fig.Q, a dip 
structure in the electronic density at a more negative gate 
voltage and corresponding current peak (dotted (green) 
and dashed-dotted (blue) lines in Panel (a-b) of Fig. ([9])) 
is visible. Actually, the above structures are broadened 
due the reduction of the quality factor as a function of the 



magnetic field. In the limit of very large magnetic fields, 
if we keep fixed the amplitude of the external antenna, 
the fine structures outlined above are completely washed 
out due to the decrease of the device quality factors. 



V. CONCLUSIONS AND DISCUSSION 

In conclusion, we have studied a CNT-based electronic 
transistor in the presence of an external magnetic field 
perpendicular to the current flux. We were able to show 
that the application of a transverse magnetic fleld modi- 
fies the bending mode CNT dynamics giving an enhanced 
damping as well as a noise term originating from the elec- 
tronic phase fluctuations induced by the displacements of 
CNT itself. 

The effective force acting on the resonator is modified 
by a pure nonequilibrium correction term proportional 
to the magnetic field as well as to the electronic current 
fiowing through the CNT. Damping and diffusive terms 
are both modified by quantum electronic current-current 
as well as density-current fluctuations corrections whose 
strengths are quadratic and linear in the magnetic field, 
respectively. 

Within our model, a quadratic dependence of the de- 
vice quality factor Q on external magnetic field strength, 
experimentally observed in Ref. [1 , naturally emerges. 
This behavior is understood in terms of a back-action 
of quantum electronic current flow fluctuations on the 
bending mode dynamics. A systematic study of device 
quality factor as a function of gate and bias voltage in the 
presence of the magnetic field has also been performed. 
All results are discussed observing the average charge and 
electronic current variations with respect to gate volt- 
age applied to the device and can be summarized as fol- 
lows. At a fixed electronic conducting regime, if negative 
charge variations and positive current variations occur, 
one has an enhanced damping reducing the quality fac- 
tor of the device. Vice-versa, negative charge variations 
and negative current variations reduce damping with a 
consequent increase of quality factors. 

We also show that, when the device is driven far from 
equilibrium, one can tune CNT-resonator frequencies by 
varying the external magnetic field: the peculiar (single 
or double dip) features in the CNT-resonator resonance 
frequency, obtained in different conducting regimes for 
the device, get distorted and acquire, in the limit of large 
magnetic field, a peculiar dip-peak structure that could 
be experimentally observed. 

Finally, when the device is actuated by an external 
antenna at fixed frequency and amplitude, the device 
current-gate voltage response is modified by fine struc- 
ture features any time the mechanical resonance with 
the proper nanotube oscillation frequency occurs. These 
structures can be tuned as a function of the external field 
and could be experimentally observed. In this sense, we 
have shown that, only exciting the CNT motion with 
application of an external radio-frequency antenna, one 
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can observe a magnetic field dependence of the electronic 
current. 

We point out that throughout this paper we do not 
take into account of a magnetic field with a component 
longitudinal to the CNT-resonator. This issue has been 
recently addressed in Ref.[14j and explained in terms of 
a more sophisticated theoretical schematization of the 
CNT-resonator electronic structure which has a cylindri- 
cal quasi-one dimensional shape. 

We end this section noting that it could be of outstand- 
ing interest to study the possibility to include quantum 
corrections to the oscillator dynamics as well as spin de- 
grees of freedorrPIl and electron-electron interaction ef- 
fects in the low bias regime. In particular, it has been 
recently proposed to study the CNT bending mode dy- 
namics by employing the spin-orbit couplin g b etween a 
single spin and nanomechanical displacemenlP^lHSl in the 
presence of a magnetic field. Quantum corrections be- 
comes important when the resonator and electronic time 
scales are of the same order on magnitude. In this direc- 
tion, it was shown in Ref . |30j that a magnetic field ap- 
plied perpendicular to the CNT results in negative mag- 
netoconductance due to quantum vibrations of the tube 
inducing an Aharonov-Bohm-like effeclP^on the electrons 
crossing the device. Work in this direction is in progress. 
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Appendix A: Current-current and density-current 
fluctuations 

In this appendix, we illustrate how the calculation of 
the force-force fiuctuation (Eq.(38)) can be performed 
with the nonequilibrium Green function approach. In 
particular, here we show how one can calculate in the 
adiabatic approximation the current-current 



S{t,t') ^ {SI{t)SIit')) (Al) 

and the density-current 

M{t, t') ^ {[5n{t)6i{t') + 5i{t)5n{t')]) (A2) 

fiuctuation terms appearing in Eq.(|38|. Recalling that 
/ ~ {II — Ir)/'2,; it is easy to see that S{t,t') is made of 
three contributions 



S{t,t') = ^[SL{t,t') + SR{t,t') + SLB.{t,t')] 



where 



SL{t,t')^{dlL{t)6lL{t')), 

SR{t,t')^{siRit)siR{t')), 
SLR(.t,t')^-{{siL{t),siR{t')}). 



(A3) 

(A4) 
(A5) 
(A6) 



The density-current fluctuation term M(t,t') is given by 
Mit,t') = ^[Mn{t,t') - MLit,t% (A7) 



ML{t,t') = {{Sh{t),Sh{t')}), 
Mjiit,t') = {{Sin{t),Shit')}). 



where 

(A8) 
(A9) 
(AlO) 

Above, {A, B} = AB + BA is an anti-commutator. In 
this appendix we limit to calculate S'L(i, t') and M]^{t, t') 
in the adiabatic approximation, since for the other fiuc- 
tuation terms the derivation is similar. 

We recall the expression for the current operator 
(through the left barrier 

iL^'^Y.^VLAd-Vl^.d^c,], (All) 

k 

We define SlL{t) — Ibit) — (II), and plan to evaluate the 
correlation function (we set V^^k = V^) 

SLit,t') = l{{siLit),siL{t')}) 



{{iL{t),iL{t')}) 



1 /ie\ 2 



(?)■ 



fe.fe' 



VkV,,{cl{t)d{t)ci,{t')d{t')) 



-Vkv:,{cimt)d\t')ck>it'))^ 
-v:vk^{d\t)ckit)cUt')d{t'))- 

+ v:v,*,{d\t)ckit)d\t')c,,it')) 



h.c 



(A12) 

The Fourier transform of S is called the noise spectrum; 
in what follows we shall be particularly concerned with its 
zero- frequency component, S{0) = J d{t—t')S{t—t') that 
is the relevant quantity in the adiabatic expansion. In 
order to evaluate the (nonequilibrium) expectation val- 
ues occurring in Eq.(Al2) in a systematic way, we first 



define the following contour-ordered two-particle Green 
functions (we follow Ref. [47 ) 

Gfir,r') = ^'{Tccl{r)dir)cl,iT')dir')) 

Gf{T,T') = z2(rc4(r)d(r)dt(r')cfe,(r')) 

Gf{r,r')^t^{TcdHr)c,{r)cl{T')d{r')) 

Gt{T,T') = l^{Tcd\T)ck{T)d\T')ck'{T')) (A13) 

The nonequilibrium noise correlator is then given by 



k,k' 



1414'Gf >(i,t') + 



VuV:,Gt>{t.t') ~ v:Vk,Gf-''{t,t') + 



cd.> , 



+ V^v:.Gt>{t,t') 



h.c. - (It 



(A14) 
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where G'f'^ {t,t') are the greater than components of 
the contour-ordered counterparts G'('^{t,t') defined in 
Eq.(Al3|. In the adiabatic approximation, we consider 



the zero-order terms of ah Green's functions G. After 
lengthly but straightforward calculations, starting from 
Eq.(|Al4| we get (we follow Ref. 47J) 



S{t, t') ^ [SL{t, t') + SR{t, t') + SLR{t, t')]/2 

^ D"{x)S{t-t'), (A15) 



where (x) is given by Eq.(43) of the main text 



Eq.(|43| is a well-known, and important result. The first 
term accounts for thermal noise (i.e., it vanishes at zero 
temperature), while the second term is a nonequilibrium 
term (shot noise), which vanishes at zero bias. 



For the mixed current-density contribution in the sec- 



ond line of Eq.(38), we have (for the left lead) 



ML{t,t') - {{SlL{t),Sn{t')}) = 
^{{lUt),n{t')})-2{hr{nr 



le 

h 



E 



v,Aclmt)d\t')d{t')) 



-v:jdHtMt)d\t')dit')) 



h.c.~2(ir)HhY. 



(A16) 



As for the the current noise spectrum S", in what fol- 
lows we shall be particularly concerned with the zero- 
frequency component of Mi(t,t'), Ml(w = 0) = J d{t — 
t')ML{t — t') that is the relevant quantity in the adiabatic 
expansion. In order to evaluate th e (non equilibrium) ex- 
pectation values occurring in Eq.(A16l in a systematic 



way, we first define the following contour-ordered two- 
particle Green functions 



G^Ti^^^') = ^'{Tccl{r)d{T)dHT')dir')), 

GII\t,t') = ^\Tcd\T)c^,(T)d\T')d{T')). 



(A17) 



In terms of the previous Greens function in Eq.(A17l, 



is then given by 

^^^(M') = f E 



(t,t') 



V,,G','f^>M + 
+ h.c.-2{lLf{nf, 



(A18) 



where Gf^^"^'^ {t,t') are the greater than components of 
the contour-ordered counterparts G^£'^{t,t') defined in 
Eq.(Al7|. Following the same reasoning as previous sec- 



tion one can show that 



dTlG{t,Tl)^L{Tl,t') 



G<{t',t) 

+G<{t,t' 

G>{t\t) 



f driSi(t,Ti)G(Ti,t') 
Jc 

dTiG{tM)^L{Tut') 



[ dTl^L{t.Tl)G{T^,t') 

Jc 



(A19) 



where is the self-energy contribution due to the cou- 
pling to the left lead and the integration is extended along 
the Keldysh contour. The function, 



[ dTlG{t,Tl)EL{Tl,t') 

Jc 



(A20) 



can be calculated using Langreth's rule^^, giving 

/(i,<')= / dhG'-{t,h)^<{t,,t') 



+ G<(t,ti)S2(ti,0, 



(A21) 



where is the advanced component of the left lead self- 
energy. In the adiabatic approximation, we consider the 
zero-order terms of all functions G and S. After lengthly 



but straightforward calculations, starting from Eq.(A19) 
we get 

M{t,t') = [MR{t,t') - ML{t,t')]/2 = D"^{x)5{t-t'), 

(A22) 



The nonequilibrium current-density noise correlator M where D^^{x) is given by Eq.(42|. 
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